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We analyze populations of Kuramoto oscillators with a particular distribution of natural fre- 
quencies. Inspired by networks where there are two groups of nodes with opposite behaviors, as 
for instance in power-grids where energy is either generated or consumed at different locations, we 
assume that the frequencies can take only two different values. Correlations between the value of 
the frequency of a given node and its topological localization are considered in both regular and 
random topologies. Synchronization is enhanced when nodes are surrounded by nodes of the op- 
posite frequency. The new theoretical result presented in this paper is an analytical estimation for 
the minimum value of the coupling strength between oscillators that guarantees the achievement of 
a globally synchronized state. This analytical estimation, which is in a very good agreement with 
numerical simulations, provides a better understanding of the effect of topological localization of 
natural frequencies on synchronization dynamics. 



PACS numbers: 05.45.Xt, 



i.75.Fb 



> 
00 

in 
o 

o 



X 



I. INTRODUCTION 



During the last decade, the realization that many ob- 
served complex systems in nature and society are intrin- 
sically related to nontrivial patterns of interactions has 
created a great interest in what we nowadays know as 
complex networks theory P, 0, d, d, Once the ba- 
sic statistical properties of the individual constituents 
(degree distributions, clustering coefficients, distances, 
diameters and so on) were analyzed in detail (see Ref. 
for a survey of measurements) the attention turned 
to the correlations among neighbors. Examples of this 
are degree assortativity [7], degree-degree correlations H 
or, even, what is known as the rich-club coefficient [9|. 
From the topological point of view, at higher scales com- 
plex networks are usually organized in groups that have 
a dense internal connectivity. These groups are known as 
functional groups in biology or as network communities 
in the social sciences. The analysis of these mesoscopic 
structures has also received much attention from the sta- 
tistical physics community '}Sk illi] • It has been also clar- 
ified how topological complexity gives rise to a complex 
dynamical behavior. Complexity results not only from 
the appearance of emergent behaviors from simple dy- 
namical units, but an important part of this collective 
behavior lies precisely in the interaction patterns [1, . 
Among the many dynamical behaviors that have been 
considered in the statistical physics literature, one of the 
most intensively analyzed is that of synchronization. We 
should understand the phenomena of synchronization as 
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an emergent cooperative behavior, where units with sim- 
ilar individual behaviors develop a coordinated collective 
output, in which all units follow the same evolution in 
time 13] . One of the paradigmatic models for the expla- 
nation of this widely occurring phenomenon goes back to 
Kuramoto [ij, UBl • In his model, units are characterized 
only by their phases ipi, which evolve according to the 
equation 



- cr^aij sin((/?j - ipi). 



(1) 



Here uji is the natural frequency of unit i, a is the cou- 
pling strength, and is the connectivity matrix [aij = 1 
if i and j are connected, otherwise). This set of equa- 
tions gives rise to two qualitative behaviors that can be 
easily understood from the model. The first term makes 
the units to follow their natural frequencies oji (in general 
these frequencies are not identical, but given by a certain 
distribution). The second one makes the phases to ap- 
proach each other. If the first term dominates we have 
incoherent behavior given by a distribution of effective 
frequencies, (fi, whose precise values depend strongly on 
the network topology. 

If the second one dominates, we observe a coherent be- 
havior where all the phases approach, and the effective 
frequencies become equal, settling the system into a syn- 
chronized state. Consequently, there must be a transition 
from an incoherent to a synchronized state. This transi- 
tion is typically characterized by an order parameter r, 
defined through the equation 



(2) 



where ^I^ is a global phase (not constant) [lj,ll5|. This 
transition has been widely analyzed for symmetric distri- 
butions of frequencies for regular [l5| as well as complex 
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[iq topologies. Frequencies are usually assumed to be 
randomly distributed across the network, and no par- 
ticular consideration about the correlation between fre- 
quency and network location has been made until quite 
recently. In particular, in Ref. 17[ a uniform distribu- 



tion of natural frequencies in different topologies has been 
considered. It has been found that, when frequencies be- 
tween neighbors are negatively correlated (i.e. nodes with 
positive frequency, taken from a symmetric interval, tend 
to be surrounded by nodes with negative frequencies) 
the network is more easily synchronized. Moreover, Ref. 
[l8| studied the synchronization in Erdos-Renyi networks 
when adaptively changing the topology of the interaction 
network to enhance synchronization. One observes the 
emergence of clusters, characterized by similar values of 
nominal frequencies, which were able to synchronize in- 
dividually significantly below the global synchronization 
threshold. The resulting interaction network eventually 
changed its topology. Namely it evolved from its initial 
random structure towards a small-world pattern. 

Following this research line, we will analyze here the 
correlation (or anticorrelation) between frequencies and 
localization in the network, but in a somewhat different 
framework. We assume that the frequencies of the os- 
cillators are distributed in a bipolar way, having values 
either -1 or +1. We take -1-1 and -1 without loss of gener- 
ality. On the one hand, the fact that the distribution of 
frequencies is bipolar does not imply it to be symmetric 
in general. Nevertheless the frequency distribution can 
be shifted arbitrarily by transforming it into a comoving 
frame of reference and the interest lies only on the phase 
and frequency differences and not on absolute values. On 
the other hand, we are only interested on the relative im- 
portance of the two terms, the frequency values and the 
coupling strength; fixing one of these values (-1-1 and -1 
in our case) just fixes the time scale. This particular dis- 
tribution (which is mathematically described as the sum 
of two delta functions) can also be understood as a limit- 
ing case of a bimodal distribution of natural frequencies, 
which has been analyzed so far in populations of all-to- 
all coupled Kuramoto oscillators 0,120, HH, i-e. without 
considering network effects. 

Also from a practical point of view, this bipolar distri- 
bution can be related to the electrical power-grid, where 
energy is either generated or consumed and hence nodes 
play complementary roles. Actually, generators have nat- 
ural frequencies slightly above the network frequency 
(50 or 60 Hz) while machines have natural frequencies 
slightly below. If we consider the system in a co-rotating 
frame (at 50 Hz) it turns out that generators have posi- 
tive small frequencies and machines have negative small 
frequencies, that corresponds to the bipolar distribution 
considered in our work. As it is shown in Ref. [HI, 
Kuramoto-like equations can be derived from swing and 
power flow equations using various simplications, as for 
example, assuming uniform maschines and generators. 
Although, a new ingredient (an inertia term) enters the 
description, which can be of crucial importance (23l. [23|. 



The final goal of the powergrid is to work in a globally 
synchronized way and the lack of synchronization can re- 
sult in serious damages. 

Finally, as we will show later, a bipolar distribution 
allows one to calculate analytically the synchronization 
threshold in a regular lattice and to determine it approx- 
imately in the general case. This last point is, actually, 
a key contribution to the literature on dynamics on net- 
works, where the number of exact results is scarce. Our 
paper is organized as follows: Section |TT] introduces the 
dynamical model and the procedure used to distribute 
the natural frequencies over the network. Then, the influ- 
ence of the natural frequencies distribution on two mod- 
els of random networks (a regular random network and a 
Barabasi-Albert scale-free network) is analyzed numeri- 
cally in Sect. IIIII This study continues in Sect. IIVI bv 
carrying out analytical calculations for the case of regular 
networks and is extended in Sect. |Vl providing a good 
approximation for the case of random networks. Section 
IVII summarizes our conclusions. 



II. DYNAMICAL MODEL AND RESHUFFLING 
PROCEDURE 



We consider an undirected graph Q which is composed 
of a set of N nodes and a set of L links, and whose 
adjacency matrix is aij. At each node of the network 
there is an oscillator whose phase evolves according to 
Eq. Il]) and is, hence, coupled to a set of neighbors. We 
assign natural frequencies -1-1 to half of the population 
and -1 to the other half at random. Starting from initial 
random phases, we let the system evolve for a time long 
enough to reach a stable state in which all oscillators 
have a bounded frequency (pi (which in general is not 
constant). 

Since the in the interaction term of the oscillators equa- 
tions is an odd function and the connectivity matrix is 
symmetric, summing up the complete set of equations of 
motion we always get 



= 0. 



(3) 



Therefore, one possible solution is tpi — for all i, which 
corresponds to the synchronized state. However, in or- 
der for this state to be reached, the coupling needs to be 
strong enough to overcome the natural frequency distri- 
bution (in our case either -1-1 or -1). Hence, for weak cou- 
pling we expect the system to develop a nonhomogeneous 
distribution of effective frequencies ipi, corresponding to 
a non-synchronized state. A transition from one state to 
the other is produced at some critical value of the cou- 
pling strength, whose determination is one of the goals of 
this paper. In order to check, in a controllable way, the 
dependence of the critical value on the precise location of 
the respective oscillators, we introduce a parameter that 
quantifies this correlation. For a single node, we define 



its frequency similarity as 



where Ni is the number of neighbors of node i and Niiuji) 
specifies how many out of them have the same natural 
frequency as node i. Furthermore, we define the overaU 
frequency similarity of the network by the mean value 

Sg = {S,). (5) 

In a similar framework, Brede [l7| analyzed a network 
of Kuramoto oscillators in which the natural frequencies 
were taken from a uniform distribution in the interval 
[— 1,+1]. By switching the frequencies of neighboring 
oscillators such that they tend to be anticorrelated, the 
author showed that this procedure enhances the synchro- 
nizability of the networks in the sense that the critical 
coupling necessary for the synchronization transition is 
decreased. 

In the following we will propose a mechanism that does 
both, either maximizes or minimizes the frequency sim- 
ilarity of the network and hence allows us to study its 
impact on the network dynamics. We proceed as follows 
Ell]: 

1. First, we randomly assign the frequency +1 to one 
half of nodes and the frequency —1 to the other half. 
Then, we measure the value of network similarity. 

2. If either the current value of similarity is close 
enough to the desired value or the number of it- 
erations exceeded a predefined limit, the procedure 
is terminated. Otherwise we continue to the next 
step. 

3. We randomly choose two nodes i and j with dif- 
ferent frequencies and we evaluate whether the ex- 
change of their frequencies increases or decreases 
the value of network similarity (see Fig. [1]). When 
the new value of similarity is closer to the desired 
one, the swap operation is accepted and the algo- 
rithm continues with the processing of step 2. 



III. SIMULATION RESULTS FOR RANDOM 
NETWORKS 

Following the procedure outlined above, we have per- 
formed simulations for two types of random topolo- 
gies. On the one hand scale-free (SF) networks grown 
by the preferential attachment mechanism proposed by 
Barabasi and Albert [25j . On the other hand random reg- 
ular (RR) graphs for which the number of links of each 
node is constant but randomly distributed to the rest of 
nodes [1^ . In both cases we take N = 500 nodes and an 
average connectivity equal to 3 (for the RR graph, 3 is 
the number of links for each node). 
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Time step t Time step t + 1 



FIG. 1: Illustration of the swap operation between nodes 
i and j in the reshuffling procedure. To represent visually 
the bimodal distribution of frequencies, we assign colors to 
nodes according to their nominal frequencies Ui. Column (a) 
corresponds to the iteration step t; two nodes i and j are 
chosen at random and we compute their frequency similarity 
Si and Sj. Column (b) represent the iteration step t -f 1, 
where the frequencies of the nodes have been interchanged 
since the overall network similarity increases. 
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FIG. 2: Order parameter r (see Eq. [2} as a function of 
the coupling strength a for random allocations of natural fre- 
quencies on two types of networks. Black and red solid lines 
correspond, respectively, to a Random Regular (RR) network 
with similarity Sg = 0.44 and a Barabasi-Albert scale-free 
(SF) network with similarity Sg = 0.624 (Number of nodes 
= 500 in both cases). Each point is an average over 30 
independent simulations with initially uniformly distributed 
random phases, and the error bars show the standard devia- 
tion. 



In Fig. [2] we show the usual order parameter r as a 
function of coupling strength cr, Eq. (2), for these two 
types of networks. Unexpectedly, r values do not increase 
monotonically with a, but show a local maximum at low 
values of a (around 0.25 for the RR and 0.1 for the SF). 
This phenomenon can be explained as a result of local 
synchronization among neighboring oscillators with the 



4 



same natural frequency when <j is small. As the cou- 
plmg strength is high enough to let these oscillators to 
interact, they tend to synchronize their phases, leading 
to higher values of the r order parameter (which is a mea- 
sure of global phase similarity). This effect disappears for 
higher values of ct, as the system evolves towards a global 
synchronization. Consequently, taking into account that 
r does not clearly differentiate between local and global 
synchronization, we adopt the following order parameter: 




which is a measure of the effective frequency dispersion 
[l^]. In our case, for a coupling strength of zero the 
order parameter is 1, whereas in the synchronized state 
it is equal to 0, facilitating a clear distinction between 
these two states. 

In the three panels of Fig [31 we can see how the tran- 
sition can be very well characterized by this new order 
parameter. In the top panel we show the network with 
the random allocation of the natural frequencies, whereas 
in the lower panels we have the behaviors corresponding 
to the extreme values of network similarity. From this 
comparison, we get a two-fold message. First, SF net- 
works synchronize more easily than RR networks, which 
is in agreement with the results of Ref. [1^. Second, low 
similarity enhances synchronizability whereas high simi- 
larity makes it harder. This last conclusion is the same 
that was obtained in Ref. when dealing with global 
synchronization. The same author, however, also realized 
that high similarity can be better for local synchroniza- 
tion (29|. The phenomena of local synchronization was 
also analyzed in a temporal perspective in Ref. (30[ . In 
this paper, a population of identical oscillators gets par- 
tially synchronized at different time scales that can be 
related to topological scales. In the case presented here 
of nonidentical oscillators, this could be achieved by in- 
troducing a discrete (but large) set of natural frequencies, 
one for each cluster to be synchronized. Similar results 
(i.e., formation of cluster with similar frequencies) were 
also observed in Ref. [l^ , as the links were locally reshuf- 
fled, increasing the similarity of frequencies in the local 
neighborhood of nodes. 

Here we notice that for low similarity the two network 
types show almost the same behavior, while for high sim- 
ilarity, which is in general worse for the synchronization, 
the SF network synchronizes more easily than the RR 
one. Actually, the SF network shows only a small ef- 
fect of the frequency similarity of the network but for 
the RR network it is evident. These results can be ex- 
plained from synchronization among neighboring oscilla- 
tors. When the frequency similarity of the network is 
small the system is more easily synchronized because the 
interaction is more effective, phases of neighboring units 
can be larger. On the contrary, for a large frequency 
similarity nodes are surrounded by nodes of the same 




FIG. 3: Synchronization (averaged stationary value of the 
proposed order parameter rS) as a function of the cou- 
pling strength ct, for RR and SF networks (A'^ = 500 and 
= 3) and different network similarities, (a) original net- 
work similarity [Srr — 0.467, Ssf ~ 0.475). (b) low network 
similarity {Srr, = 0.099, Ssf = 0.144). (c) high network 
similarity (Sflij — 0.904, Ssf ~ 0.859). Each point corre- 
sponds to an average over 30 independent simulations with 
uniformly distributed random phases as initial conditions and 
the error bars stand for the standard deviation. The vertical 
dashed lines correspond to the approximate analytical solu- 
tions for the critical value as calculated in Sec. [V] 

frequency and the system is locally stable; however this 
local synchronization is responsible for the difficulty in 
achieving the global one. Oscillators only interact along 
the borders of regions of opposite frequencies and the core 
of these regions burdens the advance of the synchroniza- 
tion wave, then the larger the regions of same frequency 
the harder to be synchronized. 



IV. ANALYTICAL APPROACH FOR REGULAR 
LATTICES 

Analyzing the above mentioned effects in regular one- 
dimensional lattices has several advantages: i) we can 
construct ad-hoc networks with extreme values of the 
network similarity parameter, ii) we can perform some 
exact calculations, and iii) we can test how well the ex- 
act results serve as an approximation for general random 
networks. 

Let us first consider a 1-dimensional ring with only 
nearest neighbor connections. The two extreme cases are: 
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FIG. 4: Distributions of natural frequencies with low (left) 
and high (right) network similarity on two regular lattices 
with node degrees 2 (top) and 4 (bottom). White and black 
colors correspond, respectively, to +1 and —1 natural frequen- 
cies. See main text for details. 



find two equivalent contributions to the coupling, the two 
links connecting oscillators sitting on the border of the 
two regions (defined by the two different values of natural 
frequencies). We find that the largest phase difference 
between any pair of nodes occurs precisely at the border 
links, and then those are the corresponding terms that 
(simultaneously) will approach the minimum value of — 1. 
Then, in this case, we find that the minimum value for 
the coupling that allows synchronization is tic = -f-- 

For the sake of completeness we have also considered 
1-dimensional lattices with next-nearest neighbors: 

c) The same geometry as case a) but with 2 additional 
connections to identical oscillators for each one (see Fig. 
Sfc)). Now one has Sg = 1/2. We can also make the same 
assumption as before and the result is unchanged with 
respect to a) because the zero phase difference between 
identical units. Hence we get exactly the same critical 
value cTc = 1/2. 

d) This corresponds to add two next nearest neighbors to 

case b) (see Fig. ID^d)) and Sg = 1 — 3/N. Analogously 

to case b) we obtain summing the equations for nodes 
1 ... N. 



a) a ring with alternating natural frequencies +1 and -1 
(see Fig. |H[a)) and in this case Sg — 0. For symmetry 
reasons, we can assume that all units with natural fre- 
quency have the same frequency (y5+ and phase (^-)_ 
in the synchronized stationary state. The same holds 
for the negative ones with frequency (p^ and phase if-. 
Then it is easy to find the following relation for any node 
with positive frequency: 

(p+ = 1 + 2(7 sm{ip^ - (p+). (7) 

In the synchronized state the left-hand-side of this equa- 
tion is 0. Then the coupling has to verify 



2 sin((/j+ — if-) 

and since the sine function is bounded a > 1/2, defining 
the critical value of cr, CTc — 1/2, above which synchro- 
nization can exist (it is a necessary but not sufficient 
condition). 

b) In this case we have half the ring with positive natu- 
ral frequencies and the other half with negative natural 
frequencies, and hence Sg ^ 1 — 2/N. We number the 
positive natural frequencies in clockwise direction start- 
ing from the leftmost point (white-colored in Fig. IHb)) 
Summing up the set of equations (1) from 1 to N/2 all 
coupling terms cancel each other out apart from the first 
and the last ones: 



(pi + - ■ ■ + (pN_ = — + asm(ipt±^^ — If N_) + a sin{ip js[ — (fii). 

2 ^ 2 2 

(9) 

Again, if the system is synchronized, the left-hand side 
of the equation equals 0. On the right-hand side, we 



ipl + ■ ■ ■ + 99 JV = 

— + crsm((^«_^i - V'f -1) + 

(T sin((y5 N — ifN,) + fj sin((^iv _|_2 — ) + 

<7sm{ipN — fi) + asin{ipN-i — 'fi) + 

crsin((pAr - (10) 

As before, we can assume that the difference between 
these phases maximizes the modulus of the sine terms. 
Then in this case we have, — — N/12. 

In general, we can say that for this simple and sym- 
metric structures, the predicted values agree very well 
with those obtained in the simulations. For cases a) and 
c) the only assumption is that white and black clusters 
have common effective frequencies and phases, although 
different among them. However, in cases b) and d) we 
rely on the fact that the maximum phase difference be- 
tween any two neighboring nodes is produced at the bor- 
der between black and white clusters. This is what we 
have also observed in the simulations (see Fig. [5]). 



V. GENERALIZATION TO RANDOM 
NETWORKS 

Finally, we present a possible extension of these results 
to the general case of random networks. Let us consider 
a network with a very low similarity, so we can assume 
that a unit with natural frequency oji is surrounded by 
Ni{uj.i) > neighbors with the same frequency and by 
Ni — Ni{uji) nodes of opposite sign. We can assume as for 
case a) in the previous section, now it is an approximation 
not as exact as before, that all positive units have the 
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which contains our previous result for isolated nodes (fT2|) 
as a particular case. In this way, we can take eq. (|14p as 
a result valid for any configuration, no matter how low 
or high is the similarity. We have tested this approach 
on scale-free and regular random networks finding good 
agreement between this prediction and numerical results 
(see Fig [3]). Our result is exact only for a particular 
case of symmetric organization. However, notice that 
simple ideas providing results on regular networks are 
only rarely so easy to generalize to any kind of network 
topology. Furthermore, this analysis coincides with 
the qualitative explanation at the end of Sect. III. We 
discussed there that the global synchronization for net- 
works of high frequency similarity was burdened by the 
size of large regions of identical frequency, which explain 
why they are harder to synchronize than networks of 
lower frequency similarity. 



FIG. 5: Synchronization in 1-d regular networks (A'' — 500) 
for different degrees and network similarities. Labels corre- 
spond to those of Fig. ID Averaged frequency dispersion 
{{rui)) as a function of the coupling a. Vertical dashed lines 
indicate our analytical estimate for the minimum coupling 
strength needed for synchronization ctc (0.5 in cases (a) and 
(c), 125 in case (b) and 41.6 in case (d)). 



same frequency and ip+ and phase Lp+, and the same 
assumption for the negative ones. Under this hypothesis 
we have : 



ipi ^ uji + {Ni - Ni{u!i))(T siii{ip-^^ - ip^^) 



(11) 



and then the minimum value of the coupling to compen- 
sate the frequency \uji\ = 1 is l/{Ni — Ni{uji)). This has 
to happen for all the nodes simultaneously in the network 
and, hence, we will have to look for the global maximum: 



1 



(12) 



On the contrary, for high similarity networks we can 
proceed in a similar way as in cases b) and d) before. 
Take a group. A, of connected oscillators with the same 
natural frequency loa surrounded by oscillators with fre- 
quency —loa, and sum up all their equations in the syn- 
chronized state 



= Na ■ to A + o" 



fly sm{(pj - ipi). (13) 

ieA.j^A 



Now we can make the approximate assumption that the 
phase difference that maximizes the modulus of the sine 
terms again is at the borders between groups of different 
sign. Then the minimum value of a that can fulfill eq. 
(fT3)) is Na/La, where La is the number of links pointing 
out of the group A. Again, since this has to be verified 
for any group, the critical value of the coupling will be 



Na 

max 

A La 



(14) 



VI. CONCLUSIONS 

In this paper we have studied synchronization in the 
population of Kuramoto oscillators with a bipolar dis- 
tribution of natural frequencies. Such scenario can, for 
example, be relevant for special cases when the network 
is optimized with respect to local or global synchroniza- 
tion [17, 29j]. It can also be seen as a proxy model for 
power grid networks where there are two complementary 
populations of nodes either generating or consuming en- 
ergy, and the final goal is to achieve a complete syn- 
chronized state where electricity is distributed in a well 
defined global frequency [23 |. 

First, we have analyzed the synchronization on random 
networks under different conditions, controlling the level 
of frequency similarity. We introduce frequency similar- 
ity of a network as a measure of the frequency similar- 
ity of neighboring nodes with respect to the frequency 
they have. We find a strong impact of the topological 
distribution of frequencies on the synchronization pro- 
cess. When the frequency similarity is low, the effect of 
network topology is only marginal. When increasing the 
similarity, the role of network topology becomes stronger, 
e.g. synchronization in scale-free networks becomes eas- 
ier than in regular random networks due to the influence 
of hubs. 

Furthermore, we have gone beyond the results in (l7j 
and by calculating exactly the minimum value of cou- 
pling strength that is necessary for global synchronization 
in simple regular networks. Generalizing this particular 
case, we obtain a simple method to estimate the syn- 
chronization threshold for arbitrary network topologies. 
Altogether, it is sufficient to count the size Na of inter- 
connected clusters of nodes, which have the same natu- 
ral frequency uja and the number of links La connect- 
ing these clusters with the nodes of opposite frequency. 
Then, the maximal ratio Na/ La calculated over all clus- 
ters present in the network gives us an estimate for the 
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coupling strength sufficient for global synchronization. 
This general result computed in terms of the topologi- 
cal distribution of opposite frequencies provides an accu- 
rate estimation of a dynamical property of the complete 
network. 
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